STATISTICS OF WAVE FUNCTIONS FOR A POINT 
SCATTERER ON THE TORUS 
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Abstract. Quantum systems whose classical counterpart have ergodic 
dynamics are quantum ergodic in the sense that almost all eigenstates 
are uniformly distributed in phase space. In contrast, when the classical 
dynamics is integrable, there is concentration of eigenfunctions on in- 
variant structures in phase space. In this paper we study eigenfunction 
statistics for the Laplacian perturbed by a delta-potential (also known 
as a point scatterer) on a flat torus, a popular model used to study the 
transition between integrability and chaos in quantum mechanics. The 
eigenfunctions of this operator consist of eigenfunctions of the Laplacian 
which vanish at the scatterer, and new, or perturbed, eigenfunctions. 
We show that almost all of the perturbed eigenfunctions are uniformly 
distributed in configuration space. 



1. Introduction 

Quantum systems whose classical counterpart have ergodic dynamics sat- 
isfy Schnirelman's theorem, which asserts that almost all eigenstates are 
uniformly distributed in phase space in an appropriate sense |2H [71 [23] . In 
contrast, when the classical dynamics is integrable, there is concentration 
of eigenfunctions on invariant structures in phase space. In this paper we 
study eigenfunction statistics for an intermediate system, that of a point 
scatterer on the flat torus. 

The use of point scatterers, or (^-potentials, goes back to the Kronig- 
Penney model |13] which is an idealized solvable model used to explain 
conductivity in a solid crystal and the appearance of electronic band struc- 
ture. They have also been studied in the mathematical literature to explain 
the spurious occurrence of the Riemann zeros in a numerical experiment |8j . 
Billiards with a point scatterer have been used extensively in the quantum 
chaos literature, starting with Seba [19], to model quantum systems strongly 
perturbed in a region smaller than the wavelength of the particle. 

The flat torus is a standard example of a system for which the geodesic 
flow is completely integrable. Placing a scatterer at a point xq in the torus 
does not change the classical dynamics except for a measure zero set of 
trajectories, and gives a quantum system whose dynamics is generated by 
an operator formally written as 

(1.1) -A + a6^, 
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with being the Dirac mass at xq and a being a couphng parameter. 
Mathematically this corresponds to picking a self-adjoint extension of the 
Laplacian —A acting on functions vanishing near xq (see Section § [3] and 
Appendix §[A]). Such extensions are parameterized by a phase (j) G (— 7r,7r], 
with (p = TT corresponding to the standard Laplacian (a = in (jl.ip ). We 
denote the corresponding operator by — A^j^^^, whose domain consists of a 
suitable space of functions f{x) whose behavior near xq is given by 

(1.2) /(x) = c(cosf i^^^^ + sm^)+o(l), x ^ xo 

for some constant C. For (p = tt the eigenvalues are those of the standard 
Laplacian. For 7^ vr (a 7^ 0) the resulting spectral problem still has the 
eigenvalues from the unperturbed problem, with multiplicity decreased by 
one, as well as a new set A,^ of eigenvalues interlaced between the sequence 
of unperturbed eigenvalues, each appearing with multiplicity one, and sat- 
isfying the spectral equation 

(1-3) ^\i^n{xo)\'^{^--^- ^^^) = cotan^ 

for a certain cq > 0, where {ipnix)} form an orthonormal basis of eigen- 
functions for the unperturbed problem: —Aipn = ^ntpn- The eigenfunction 
corresponding to A G A,^ is the Green's function 

(1.4) Gx{x;xo) = {A + X)-^6x, . 

Our main result is that for almost all A S A,^, the perturbed eigenfunctions 
Ga(»; Xq) are uniformly distributed in position space. To formulate the result 
precisely, we denote by 

(15) a Jx) - ^^(^'^P) 

(1.5) 9x{x).- ii^^ii^ 

the L^-normalized Green's function: 

Theorem 1.1. Fix (j) G (— 7r,7r). There is a subset A^^oo C A^ of density 
one so that for all observables a £ C°°(T^), 

(1.6) / a{x)gx{xfdx 7^1 a{x)dx 

Jf2 area(T^) Jf2 

as A —)• cxD along the subsequence A^ 00 
Remarks: 

For the eigenfunctions of the unperturbed Laplacian, there is a variety 
of possible limits in the position representation, which were investigated by 
Jakobson [lOj. 

A result of the same nature as our Theorem 11.11 was recently obtained 
in [15] for billiards in rational polygons. There it is shown that for any or- 
thonormal basis of eigenfunctions, there is a density one subsequence which 
equidistributes in configuration space. The method of [15] adapts the proof 
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of quantum ergodicity for billiards of |24] to work in configuration space and 
inputs the theorem of Kerckhoff, Masur and Smillie [12] who showed that 
for rational polygons, the billiard flow is uniquely ergodic in almost every 
direction. Our argument here is completely different and is very specific to 
this particular model. 

A related, and in some sense complementary, issue was studied by Berko- 
laiko, Keating and Winn [2] who predict that for an irrational torus with a 
point scatterer there is a subsequence of eigenfuctions which "scar" in mo- 
mentum space, and this was proved by Keating, Marklof and Winn [11] to 
be the case assuming that the eigenvalues of the Laplacian on the unper- 
turbed irrational torus have Poisson spacing distribution, as is predicted by 
the Berry- Tabor conjecture. 

It is important to note that we (as well as [21 [11]) deal with the limit 
of large energy A — >• oo for a fixed phase 7^ tt, which is called the weak 
coupling limit in the physics literature. An interesting problem would be to 
understand the strong coupling limit, where A — >• 00 together with — t- vr 
while tan((/)/2) ^ log A, so that the RHS of the spectral equation ()1.3p blows 
up. In that range it has been argued that the spectrum displays intermediate 
statistics [20l El [3l H [JT] . 

Acknowledgments: We thank Maja Rudolph for her help with the numeri- 
cal investigation of some of these issues, and John Friedlander for discussions 
concerning sums of two squares. Z.R. was partially supported by the Israel 
Science Foundation (grant No. 1083/10). H.U. was supported by a Minerva 
Fellowship. 



2. The flat torus 

2.1. Basic setup. We consider a flat torus T-^ obtained by identifying oppo- 
site sides of a rectangle with side lengths 27r/a, 2iTa, so that = M^/27r£o 
where Cq = Z{l/a,0) © Z(0,a) is a unimodular lattice. 

An orthonormal basis of eigenfunctions for the Laplacian A on consists 
of the exponentials 

(2.1) J_e^(^'?) 

where ^ ranges over the dual, or reciprocal, lattice 

(2.2) C = {x GR"^ : (x,i) £Z, yi£Co} = {(ma,-):m,neZ} 

a 

The eigenvalues of the Laplacian on are the norms of the vectors of 
the dual lattice C Weyl's law for the torus, establishing the asymptotics of 
the counting function N{x) of eigenvalues below x, is equivalent to counting 
the number of points of the lattice C in a disk (equivalently the number of 
points of the standard lattice Z'^ in an ellipse), and therefore reads 

(2.3) N{x) = <x:^eC} = 7rx + O(x^) 
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The exponent 9 in the remainder term is known to be at least 9 > 1/4. 
The "trivial" bound on the remainder term, as the length of the boundary, 
translates into 9 < 1/2. A nontrivial bound uses Poisson summation and 
the method of stationary phase leads to 9 < 1/3. We will need a better 
bound 

(2.4) 9 < ^ 

such as the one due to van der Corput [22j . The current world record of 
9 < 131/416 + o(l) is due to Huxley [9]. 

Using the remainder term (|2.3p . we may deduce a bound for the number 
of lattice points in an annulus: Define 

(2.5) A(A,L) = G £ : A-L < < A + L} 



Then ((231) implies 

(2.6) #A(A,L) = 27rL + 0(A^) 

2.2. Multiplicities. Denote by A/" = {0 < ni < . . . } the set of norms of 
the dual lattice vectors. The multiplicities in the spectrum are 

(2.7) rein) = #{e € £ : I^P = n} 

The lattice C is rational if, after a suitable scaling, the norms are all 
rational. The norm of an arbitrary lattice vector ^ = {ma,n/a) is 



1 



(2.8) ICI^ = aW + nVa^ = ^(aW + n2^ 



a 



2 ' 



SO that the lattice is rational if and only if G Q is rational. 

In the irrational case, the multiplicities are entirely due to the reflection 
symmetries (x, y) i— )■ (ix, iby) and are given by 



(2.9) rein) 



1, n = 

2, n = a^u^ or v'^/a'^, u,v £ 
4, otherwise 



and in particular are generically equal to 4. 

For the rational case the multiplicities are complicated arithmetic func- 
tions. For instance in the case of the standard lattice £ = Z^, the multiplicity 
'"("-) = i^ii^^y) G : + = n} depends on the prime factorization of 
the integer n. In any case, it is well known that we have an upper bound on 
the multiplicities of the form (see e.g. the proof of \16\ Lemma 7.2]) 

(2.10) rc{n) < n°(^) 
The counting function of the norms is 

(2.11) Mix):=#{neM :n<x} 
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Since we have the upper bound ()2.10l) and since we know that the sum over 
n < X of the multiphcities is asymptoticahy vrx ()2.3p . we deduce a lower 
bound 

(2.12) Afix) > xi-°(^) 

2.3. Nearest neighbor gaps. If the norms are ordered by = {0 < ni < 
n2 <...}, we need to understand the spacings (or gaps) n^+i — between 
successive norms. The individual values are difficult to understand. From 
(j2.3p we certainly have n^+i — 71^ ^ n^. However one can do better by 
arguing as follows [Tj: First find the largest integer square < n^+i/a^, 
which one can do so that n^+i — a^u^ <C ^nu+i- After that find the largest 
square v'^ /a? < nk+i — a^v? , which one can do so that n^+i — a^u^ — f^/a^ ^ 
-^rifc+i — a^u^ <^ Thus we found a norm n = a^v? + v'^ /a? £ M with 

n < n^+i so that n < giving 

(2.13) nfc+i - nfc < nH'^ 

The average spacing for norms up to x is, using (|2.13p . 
/ .X 1 . X x(l + o(l)) 

and by the lower bound (12.121) we deduce that 

(2.15) 1 ^(n,-n,_i)«x°« 

Since we are dealing with averages of positive quantities, we find: 

Lemma 2.1. For almost all k, that is on a density one sequence, the spac- 
ings satisfy 

(2.16) Tifc+i - Tifc <C n°^^^ 

3. POINT SCATTERERS AND J-POTENTIALS 

In this section we review the theory of point scatterers (see [6]), with some 
details left to Appendix El 

3.1. A finite-dimensional model. We want to study the Schrodinger op- 
erator with a delta-potential on the flat Z)-dimensional torus T'^ = M.^ /2ttCq, 
where Cq C is a unimodular lattice. The operator is formally given by 

(3.1) -A + a6 

where 6 is the Dirac delta- function at the point xq. 

To make sense of the operator ()3.ip . we say that for a wave function 
■ip € C°°{T^), multiplication by 6 should give 

(3.2) 6tl; = tl;{xo)6={iP,6)6 
which we try to think of as a rank-one operator. 
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As an approximation, it is useful to examine a finite-dimensional model: 
a rank one perturbation of a self-adjoint operator Hq on a finite-dimensional 
Hilbert space Ti of the form 

(3.3) H = Hq + aP^, PyU = {u, v)v, u G Ti 

where ^ v € T-L and a 7^ is real. Let (pn be an orthonormal basis of 
Ti consisting of eigenvectors Hq , with eigenvalues e^: HQCpn = en4'n- The 
eigenvectors (pn of Hq with {(j)n,v) = are clearly still eigenvectors of H. 
The new part of the spectrum is given by: 

Lemma 3.1. The eigenvalues E ^ Spec {Hq) of the perturbed operator H 
are the solutions of the equation 

(3.4) (^^E-Hq)-\,v) = - 

a 

or equivalently 

E%^ = - 

^ ' ^ E-En a 

n " 

with corresponding eigenfunction 

(3.6) u={E-Hq)-^v 

Proof. We rewrite the eigenvalue equation Hu = Eu for H in the form 

(3.7) {E — Hq)u = a{u,v)v 

li E ^ Spec(-ffo) then necessarily {u, v) ^ and we find that 

(3.8) u = a{u, v) {E - Hq)-^v 
Thus up to a scalar multiple 

(3.9) u = {E- Hq)~^v 
Substituting (|3.9p in the eigenvalue equation (|3.7p gives 

(3.10) V = a{{E - Hq)'^v,v)v 
that is 

(3.11) {{E-Hq)-\,v) = - 

a 

Expanding v = J2ni''^ ^ ^n) 4'n in terms of the normalized eigenvectors cpn of 
Hq gives 

K^,<^n)P 1 



(3.12) 



E - en a 



Conversely, if E ^ Spec(i?o) and (j3.4p holds, take u = {E — Hq) as in 
^Ml- Then 

{H - E)u = {Hq - E)u + a{u, v)v 

Since {Hq — E)u = {Hq — E){E — Hq)~^v = —v and a{u,v)v = a{{E — 
Hq)-^v, v)v = vhy ([13]) we find that Hu = Eu. □ 
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Now take for Hq the free Schrodinger operator Hq = —A acting on 
C°°(T^) C L2(T^), and v = 5^^. Then the eigenfunctions with eigen- 
values E ^ Spec {Hq) are the Green's function 

(3.13) Gij(.;xo) = (S-Fo)~i<5,o 

However the sum ()3.5p diverges; indeed, in that case the RHS of ()3.5p equals 

|^n(xo)P 



(3.14) 



E 



which is divergent in dimension D > 1, by Weyl's law. In fact (j3.14p is 
just the Green's function evaluated on the diagonal, which is divergent in 
dimension D > \. Thus one needs a regularization procedure. 

3.2. Regularization. One regularization procedure is through the theory 
of self-adjoint extensions. A succinct account of this procedure is given in 
[B]. For the reader's convenience this will be reviewed in Appendix § [Al 
One starts with the standard Laplacian A, for which an orthogonal basis 
of eigenfunctions are the exponentials e*^^'^\ ^ in the dual lattice C with 
corresponding eigenvalue \^'\^ ■ The idea is that for functions vanishing at 
the point xq all candidates have to coincide with the unperturbed operator 
Hq = —A. We want to extend it to a bigger space. Denoting by —/^xo 
the unperturbed operator restricted to C^(T-^\{xo}), one finds that the 
adjoint has as its domain Dom(— A*|j) the Sobolev space H'^{T^\{xq}). In 
dimension D = 2, 3 this equal^ the space of / G L^(T-^) for which 3 A G C 
s.t. 

(3.15) A/-^4o G^'(T^) 

and for such /, there is some S G C so that for x near xq, 

(3.16) /(x) = ^G(^)(|x - xol) +B + o(l), 
where 

(3.17) c^mJa;-^'-- "=1. 

One finds that in dimension D = 2, 3 there is a one-parameter family of 
extensions, parameterized by (/) G M/27rZ ~ U{1), denoted by —A^^xq with 
domain given by / G Dom(— A*^) for which there is some a G C with 

(3.18) /(x) = a ^cos (|x - Xol) + sin^^ + o(l), X — )• xq 

The action of — A^ on / satisfying (|3.15p is then given by 

(3.19) - A4>,xJ = -A/ + A6x, 



^In dimension D — 1 one wants A/ — co6 — ci&' G L^(T^), while in dimensions D > 4, 
Hq is essentially self adjoint and there are no self-adjoint extensions. 
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The divergent equation ()3.5p is replaced by the convergent equation for the 
new eigenvalues A ^ cr{A): 

where 



We can rewrite ()3.20p as 



ie|2-A lei^ + ij 2 

1 

(3.21) ^^^(n)/^-^^| =cotan^ 

I n — A n"^ + 1 I 2 

where 

rc{n) = j^{n=\i\^ -.i^C) 

and A/" is the set of norms of vectors in L (without multiplicity). The 
corresponding eigenfunction is a multiple of G\{x\xq) = (A+A)~^(53;(,, which 
has the L^-expansion 

(3.22) G>(...o) = -^E °"'f;p'"r°" . 

As may be seen from (j3.2ip . the new eigenvalues interlace between the 
sequence M = \n\ < 712 < . . . } of norms, see Figure [TJ We may thus label 
the perturbed eigenvalues A^ = A^ so that 

(3.23) ni < Ai < n2 < A2 < • • • < rifc < A/c < nfc_|_i < . . . 



3.3. The density of states. The density of states of the perturbed eigen- 
values depends strongly on the particular torus, that is on the lattice it 
determines. 

For the standard lattice Z^, a theorem of Landau |14j asserts that 

(3.24) #{n G AA, n < x} ~ B J^ , 

ylogx 

where B = -^n(-'- ~ p^^)~^^^ = 0.764. . . , the product over primes p = 
3 mod 4. Consequently we deduce a form of Weyl's law for the perturbed 
spectrum of A^^j.^ for the standard lattice: 

X 



(3.25) #{A G A^ : A < x} ~ B- . 

Vlogx 

In the irrational case, the multiplicities are typically 4, see (12. 9p . Then 
Weyl's law in those cases would read as 

(3.26) #{A £A^:X<x}^jX 
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Figure 1. A plot of the spectral function on the LHS of the 
eigenvalue equation (|3.2ip for the standard lattice C = 1? . 
The intersections of the plot with the x-axis are the perturbed 
eigenvalues corresponding to </> = 0, alternating with the 
norms. 



4. The norm of G\ 

We will need a lower bound on the L^-norm of the Green's function Gx. 
We are able to get a good bound for a sub-sequence of density one. To define 
this subsequence, we recall our discussion of the gaps between consecutive 
norms. 

According to Lemma l2.H for almost all k we have 

(4.1) nfc+i - nfc <C n°^^^ 

We define the set A^ C A of eigenvalues (using the labeling (I3.23P ) so 
that (gH) holds: 

(4.2) Kg = {\k e A : nfc+i - < ^} 

By the discussion above, this is a set of density one in A (and conjecturally 
all of A). 

Lemma 4.1. For A G A^ (i.e. for almost all \), we have 

1 



(4.3) I|Ga||»^,(i) 

Proof. Let = nk{X), nfc+i = nfc+i(A) be consecutive norms so that < 
A < ri/k+i. Then trivially 
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Since for A G Ag we know that Uk+i — <^ n 



^ follows. □ 
It is natural to conjecture that (|4.3p holds for all A. 

5. Truncation 

For L > let ^(A, L) be the set of lattice points in the annulus X — L < 
|xp < A + L: 

(5.1) A{X,L) = {C€C:\\C\^-X\<L} 

We denote the truncated Green's function by 



(5.2) 



— y 



47r2 



exp(i^ • {x - xq)) 



?eA(A,L) 

Let g\ and g\ l be the L^-normalized Green's function and its truncation: 



(5-3) ffA = TT7TTr> 9x,L II 

II^aN II^a,l|| 

We have the following approximation. 
Lemma 5.1. Let L = X^ , 6/2 < 5 <1. As A — )■ oo along Kg, 

Wqx - 9\,l\\2 0. 

Proof. Note that 



G 



\G 



A|l2 



\G 



A,L|l2 



< 



\G\ - Gx^lW^ 



\G 



A|l2 



(5.4) 

We have 
(5.5) 



+ I|G'a,l||2 

< 2 



1 



G'a||2 

Gx — Ga,l||2 



\G 



A,L|l2 



\Gx - G'a,l||2 



1 



IIGaII, 



167r4 ^ ^ (|^|2 - A)2 

II5P-A|>A« ^ 



We recall how to evaluate lattice sums using summation by parts: 
Let ni < n2 < • • • be the set of norms, and 

(5.6) N{t) = rcink) 

nk<t 

Then for a smooth function /(t) on M we have 
(5.7) 

/del') = NMfins) - N{nA)f{nA+i) - / nt)N{t)dt 
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Now use the lattice count with remainder (12 .Sp to get 

(5.8) /(I^l')=^/ /W^^ 

+ 0{n%f{nB) + n'^^JiuA)) + 0{ / \nt)\t'dt) 

Applying (jS.Sp with f{t) = — A)^, once with ua = ni and ub < 

X — L < UB+i and then with ha-i < X + L < ua < nA+i and = oo gives 

(5.9) ||G'a-Ga,l||2« ^ + ^ 



Since for A G we have UGaIP ^ 1/X°^^^ by Lemma |4.H we find 

^X" 



which tends to zero if 5 > 0/2. □ 

Consequently we may study the matrix elements by replacing gx by the 
truncated version g\^L- 

Lemma 5.2. Let f e C°°(T2) and L = X\ 9/2 <6 < 1. We have 
I if 9x,9x) - if 9x,L,gx,L) I 

as A —)• oo along Ag. 

Proof. Let / G C°°(T2). We define the multiplication operator M/ : L'^{T'^) - 
L2(T2) by 

Mf{g) = fg. 

Since Mj is a continuous operator on L^(T'^), we have that 

lls-A - 9x,l\\2 

for A G Ag implies 

\\Mfigx- gx,L)\\2^0 

and hence 

I {Mfgx,gx - g^Li I < ll^/IL II^a - 9x,l\\ 0. 
If we repeat this, where we switch gx and gx,L^ we obtain 
I {Mfgx,9x) - {Mfgx,L,9x,L) I ^ 0. 

□ 
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6. Exceptional eigenvalues and a Diophantine inequality 

Fix a nonzero vector {) ^ Q & C, and 5 E (|, ^ — 0) (such 5 exists because 
9 < 1/3). Let Si^ be the set of vectors satisfying 

(6.1) S^ = {v^C:\{r^,C)\<\r]n 
We define a subset C A of eigenvalues 

(6.2) A^ = {Ag A:^(A,A^)n5^ = 0} 

(recall that A{X,L) are the lattice points rj £ C in the annulus (jS.ip ). Our 
goal in this section is to show that 

Proposition 6.1. A^ is a subset of density one in A. 

Proof. Let 

(6.3) B(; = A\A^ = {A e A : A{X, A^) n 5^ ^ 0} 

We will show that has density zero in the set A of all perturbed eigen- 
values, in fact 

(6.4) #{A G : A < X} < 



ICI 



for 6' = 1/2- 9-6 >0. 
We first show that 

Lemma 6.2. 

(6.5) E : < X} « 



ICI 



Proof. Introduce cartesian coordinates with one of the axes in the direction 
of the vector by writing every x E as 

C 

(6.6) X = u— + V- 



ICI IC^I 

where if C = (oj*) f ) then (^-^ = {—^,pa) is a vector orthogonal to C- In these 
coordinates, 

(6.7) {x,C)=u\C\, |2;|2 = m2 + ^;2 

and our set of lattice points is thus contained in the rectangle 

(6.8) n = {u-^^+v-^^:\u\<^, |H<Xi/2} 

Now estimating the number of lattice points in a rectangle is a simple matter: 
Putting a fundamental domain = {{ax,y/a) : < x,y < 1} for the lattice 
(which has unit area for the case at hand) centered around each lattice point 
in TZ, we get a figure whose area is the number of lattice points in question, 
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and which is contained in a shghtly bigger rectangle whose dimensions are 
expanded by the diameter dc = yaF+l/a? of F: 

(6.9) = {u-^^+v^^:\u\<^^+dc, \v\<X^/^ + dc} 
Thus we see that #£ n 7^ is bounded by the area of 7^+, which is: 

vS ^1^1/2+5 

(6.10) area 71+ = 2(— + dc) x 2{X^I^ + dc) = + Oc{X^'^) 
Therefore 

(6.11) #{r/ G 5c : |r?|2 < X} < area7^+ = + Oc{X^'^) 

as claimed. □ 

Next we define J\f(^ C A/" to be the set of norms |r/p of ?7 G 5^^, without 
multiplicities. We clearly have 

TAl/2+5 

(6.12) #{n G AA^ : n < X} < G : |r?|2 < X} « 

by Lemma |6.2[ 
We have a map 

(6.13) l:Bc_^Mc_ 

defined by i{X) being the closest n G AA^ to A; if there are two such elements, 
i.e. n_ < A < n+ with n± G and n+ — A = A — n_, then set t(A) = n+. 
Thus we get a well defined map, whose fibers satisfy 

r^{n) C {A G A : 37/ G SqC] A{\\^),\y]\'^ = n] C A n [n - 2n^ n + 2n'^] 

for n ^ 1. 

Since A is interlaced between the norms A/", we have 

(6.14) #An [n-2n^n + 2n^] < ^ r£(m) = #yl(n, Sn"^) 

n— 3n'^ <m<n+3n* 

which is the number of lattice points in an annulus. By (|2.6p . 

(6.15) ij^A{n,2,n^) <^n^ + 
and hence (since 5 < 9) 

(6.16) #L-^{n)^n^ 
We thus find 

#{A G i?c : A < X} = ^ #6-i(n) 

(6.17) 

T^l/2+5+6» 

< X''#{n G AAf : n < X} < ^ 
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That is 

(6.18) #{A e : A < X} « 

with 5' = 1/2 - - 5 > 0. □ 

7. Proof of Theorem 11.11 

7.1. Fixed observables. Fix a nonzero vector C, ^ C and recall the defi- 
nition (j6.2p . (j4.2p of the sets of eigenvalues and A^; both are of density 
one in A and thus 

(7.1) A,,^:=A,nA^ 

is still a set of density one in A. 

We will show that {e^gx^gx) — )• as A — t- cx) along A^ ^^. By Lemma |5. II it 
suffices to show: 

Proposition 7.1. Take A G A^^^ and L = , 6 G (|,^ — 0). Fix nonzero 
C e C. As X^O while X G A^^^, 

(7.2) {e<;gx,L,gx,L) ^0 
Proof. For ^ G £ and A G Ag^(^ define 

The L^-norm of the truncated Green's function Gx^l {L = X^) is given by 

ieA{x,L) 

and hence 

(7.4) (^e^Gx,L,Gx,L) = Yl^ ^ <0c{^ - ■ 

ieA{x,L) 

Cauchy-Schwarz gives 

\{ecGx,L,Gx,L)\^ <\\Gx,l\\1 E ^(^ " O'- 

5eA(A,L) 

Note that for A G A^^^, 

Ke-C)l«^. 

Indeed 

(7.5) |^_^|2_A = |e|2-A-2(e,C) + |Cl' 

and since for A G A^ and ^ G A{X,L) we have |(C,C)I > ICP^ ~ L we find 
that 

(7.6) lie - CI' - A| > 2L(1 + 0(1)) -L- Id' » L 
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Using ()2.6p gives the bound 



(7,7) c((-Cf«l^^«-^ 



(recall 9 > 1/4) so that we find 

(7.8) {ecGx,L,Gx,L) <.\\Gx\\^ 

The lower bound ^ 1/A°(^) of Lemma |4. II implies 

for A € Ap^(;, which tends to zero since 6 > 6/2. □ 

7.2. A diagonalization argument. We have shown that for each ^ C G 
£, there is a density one subset Ag^,^ of eigenvalues so that {e,^g\,gx) — )• 
as A — )• oo along Ag^^. It remains to see that there is a density one subset 
Aoo C A so that for every observable a G C°*^(T^), we have 

(7.9) {ag\,gx)-^ 7^/ a{x)dx 

area(T^) Jj2 

as A — )• oo along Aoo- We recall the argument, which can be found e.g. in [7J. 
For J > 1, let Aj C A be of density one so that for all |C| < ^, {^c^Qx-, dx) 
as A — )> oo along A^^^, and in particular for every trigonometric polynomial 
Pj{^) = T^lciKjPC^d^) we have 

(7.10) {Pjgx,gx)^ ^/ PAx)dx 

area(Jl^) Jj2 

We may assume that Aj_|_i C Aj for each J. Now choose Mj so that for all 
X > Mj, 

and let Aoo be such that Aoo n [Mj, Mj+i] = AjH [Mj, Mj+i] for all J. Then 
Aoon[0, Mj+i] contains Ajn[0, Mj+i] and therefore Aoo has density one in A 
and ()7.10p holds for A G Aoo- Since the trigonometric polynomials are dense 
in C°°(T^) in the uniform norm and the probability measures \g\{x)\'^dx are 
continuous with respect to this norm, we find that (j7.9p holds. 

Appendix A. A rigorous description of the point scatterer 

Denote by 6xq the Dirac distribution at xo € T^. We are interested in 
solutions to the equation 

(A.l) (A + A)/ = 4o, f(^C^iT^\{xo}) 
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and its association with the eigenfunctions and eigenvalues of a family of 
self-adjoint operators. Consider the domain of C°°-functions which vanish 
in a neighborhood of xq 

Do = C^{T^\{xo}) 

and denote —^^xo = —^Idq- The operators associated with equation ()A.ip 
form the family of self-adjoint extensions of the positive symmetric operator 
-A^o (cf. [6], section 1, p. 277). 

We make the following conventions in the definition of the Green's func- 
tion. 

Definition 1. Denote by cr(— A) the spectrum of —A on C^(T^). We define 
the Green's function Gx{x;y) at energy A G C\cr(— A) on as the integral 
kernel of the resolvent (A + A)~^, that is 

(A.2) (A + A)-i/(y)= / Gx{x,y)fix)dx. 

In order to give a self-contained presentation of the theory of self-adjoint 
extensions, we briefly recall the standard definitions of the adjoint of an 
operator, symmetry and self-adjointness. 

Definition 2. Let H he a Hilbert space and Dom(i?) C H. Consider the 
operator B : Dom(i?) — ?• H. We define 

Dom(S*) = {yeH\3aeH: Vx G Dom(5) : {Bx, y) = {x, a)} 

then we define the adjoint B* of B as the map B* : Dom(i?*) — )• H, 

B*y = a. 

We call B symmetric if 

\/x, y G Dom(S) : (Sx, y) = {x, By) . 
We call B self-adjoint if B is symmetric and 

Dom(S) = Bom{B*). 

We have the following well-known results from self-adjoint extension the- 
ory which we summarize briefly. Proofs can be found in lLB\, Chapter X.l. 
We will give the relevant references for each lemma. 

Definition 3. Let i? be a densely defined symmetric operator on a Hilbert 
space H. Denote its adjoint by B* . Let r/ G C \ M. The deficiency spaces of 
B at rj and fj are defined as 

(A.3) DriiB) = kei{B* - rj}, Dr,{B) = ker{B* - fj}. 

We refer to the members of a basis of a deficiency space as deficiency ele- 
ments. 

The following lemma is proven as part of Theorem X.l, p. 136 in |18j . 

Lemma A.l. As a function of rj, dim 1)^(5) is constant on the upper 
(lower) complex half-plane. 
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We proceed with the definition of the deficiency indices of a closed sym- 
metric operator which indicate if the operator can be extended to a self- 
adjoint operator. 

Definition 4. If dimD^(i?) = m and dim Dfi{B) = n for nonnegative 
integers m, n, we say that the operator B has deficiency indices {m,n). 

The following lemma is Corollary (a), p. 141 in |18] . 

Lemma A. 2. B has deficiency indices (0, 0) if, and only if, B is self-adjoint. 

If the deficiency indices are nonzero and equal, then a family of self- adjoint 
extensions exists and can be constructed as follows. This lemma combines 
Theorem X.2, p. 140 and Corollary (b), p. 141 in [18] . 

Lemma A. 3. If a closed positive symmetric operator B has deficiency in- 
dices {n,n), n>l, then for each unitary map U : D\{B) — )• D-\{B) there is 
a self-adjoint extension Bjj : Du — )• H , where 

(A.4) Du = {f = g + h + Uh\{g,h)e Dom(S) x D;{B)] 

and Bjj acts as follows 

(A.5) Buf = Bg + ih-iUh. 

The operator Bjj has deficiency indices (0,0). Conversely, every self-adjoint 
extension of B is of this form. 

We apply Lemma lA. 31 to construct a one parameter family of self-adjoint 
extensions of the operator —Ax^. Denote the domain of the closure of — Aa,^ 
hyDo- 

Lemma A.4. The operator — A^^g has deficiency indices (1,1). The corre- 
sponding deficiency elements are the Green's functions Gi{x,xo), G-i{x,xo). 
The self-adjoint extensions of —A^q are given by the one parameter family 

(A.6) -A^:D^^ L^{T^), ^ € (-^, vr] 

where 

(A.7) D^ = {g + cG; + ce''^G_i : 5 G i^o, c G C} 

and 

(A.8) - A^f = -Ag + ciGi - ce^^iG_i. 

Proof. By definition of the Green's function we have that 
(A.9) ker{(A,J*±i} = £{G±i}. 

Hence —A^^ has deficiency indices (1,1) and we may apply Lemma lA. 31 to 
obtain the result. □ 

Remark A.5. Note that —A^^ recovers the Laplacian on C°^(T^). 

Next we derive an equation for eigenvalues of the operator — A^, (p G 
(-7r,7r). 
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Lemma A. 6. Let ip £ (— 7r,7r). We have that A ^ a{—A) is an eigenvalue 
of —Aip if, and only if, 

(A.IO) (^2^ - ^T^) = ^0 



where 



°° S|f|4 + 1- 



T/ze corresponding eigenfunction is a multiple of Gx{x;xq). 
Proof. Let / G and ||/||2 = 1. Then / must be of the form 
(A.ll) / = 5 + cGi + ce'^G_i, g€Do,ceC. 

Let us first assume that A ^ (t(— A) is an eigenvalue of — A^. We have 
(A.12) = (A^ + A)/ = (A + A)(7 + c(A - i)Gi + ce''^(A + i)G_i. 
We apply the resolvent (A + A)^-*^ to both sides to obtain 

(A.13) o = . + .f^O, + »l-^0_,. 

In view of the iterated resolvent identity 
(A.14) (A^i) ^ ^ ^ ^ 



A + AA±i Azbi A + A 
we can rewrite equation ()A.13P as 

(A.15) = g{x) + c(Gi - Gx){x, xq) + ce*'^(G_i - Ga)(x, xq). 

In particular 

(A.16) f = g + c(Gi + e*<^G_i) = c(l + e''l')Gx 

and so / is a multiple of Gx . 

If we now take the limit x — t- xq on the r.h.s. of (lA.lSh we obtain 

(A.17) = lim (Gi - Ga)(x, xq) + e'^ lim (G_i - Ga)(x, xq) 

and note that A ^ cr(— A) implies c 7^ so we may drop the constant. A 
simple rearrangement of this equation yields 

(A.18) tan((^/2) lim Im Gi(x, xq) = lim (Ga - Re Gi)(x, xq). 

X^Xq X^Xq 

In order to obtain equation (jA.lOp we require the following L^-identity for 
the Green's function Ga on 

(A.19) G,(.,.„)= 1 J^2E(|-(^L_M, 
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We rewrite the r.h.s. of ()A.18p as 

(A.20) = - 4^ E - ^o)) {^2^ - ^1^} 

-y 

L7r2 



4vr2|l^l|e|2-A |e|4 + l 
Finally, note that 

Jim^ImGKx,xo) = -^|:^^. 

To see the reverse implication assume that A solves equation (jA.lTp . a 
rearrangement of equation (lA.lOp . The r.h.s. of equation (]A.17p has singu- 
larities at points which are in (t(— A), hence A ^ (t(— A). We define 

h{x) = (Ga - Gi)(x,xo) + e''^(GA - G_i)(x,xo) 
and observe that 

(A.21) (1 + e''^)GA = /a + Gi + e''^G_i G 

because equation (IA.17P implies fx{xQ) = 0. The iterated resolvent identity 
()A.14p implies 

(A.22) (A^ + A)/a = (A + A)/a = - (A - i)Gi - (A + i)Gi 

and by the definition of the operator A^ we have 

(A.23) (1 + e'^)(A^ + A)Ga = (A + A)/a + (A - i)Gi + (A + i)G_i = 0. 
This concludes the proof. □ 



References 

[1] R. P. Bambah and S. Chowla. On numbers which can be expressed as a sum of two 

squares. Proc. Nat. Inst. Sci. India, 13, 101-103 (1947). 
[2] G. Berkolaiko, J. P. Keating and B. Winn. Intermediate wave functions statistics, 

Phys. Rev. Lett. 91, 134103 (2003). 
[3] E. Bogomolny, U. Gerland and C. Schmit. Singular statistics. Phys. Rev. E (3) 63 

(2001), no. 3, part 2, 036206. 
[4] E. Bogomolny, O. Giraud and C. Schmit. Nearest-neighbor distribution for singular 

billiards. Phys. Rev. E (3) 65 (2002), no. 5, 056214. 
[5] T. Shigehara and Taksu Cheon, Wave chaos m quantum billiards with a small but 

finite-size scatterer. Phys. Rev. E 54, 1321-1331 (1996). 
[6] Y. Colin de Verdiere. Pseudo-laplaciens 1. Annales de I'lnstitut Fourier, tome 32, no. 

3, 275-286, 1982. 

[7] Y. Colin de Verdiere Ergodicite et fonctions propres du laplacien. Comm. Math. Phys. 
102 (1985), no. 3, 497-502. 



20 



ZEEV RUDNICK AND HENRIK UEBERSCHAR 



[8] D. A. Hcjhal. Some observations concerning eigenvalues of the Laplacian and Dirichlet 
L-series. Recent progress in analytic number theory, Symp. Durham 1979, Vol. 2, 95- 
110 (1981). 

[9] M. N Huxley, Exponential sums and lattice points. III. Proc. London Math. Soc. (3) 

87 (2003), no. 3, 591-609. 

[10] D. Jakobson. Quantum limits on flat tori. Ann. of Math. (2) 145 (1997), 235-266. 

[11] J. P. Keating, J. Marklof and B. Winn. Localized eigenfunctions in Seba billiards. J. 
Math. Phys. 51 (2010), no. 6, 062101. 

[12] S. Kerckhoff, H. Masur and J. Smillie, Ergodicity of billiard flows and quadratic dif- 
ferentials. Ann. of Math. (2) 124 (1986), no. 2, 293-311. 

[13] R. de L. Kronig and W. G. Penney, Quantum Mechanics of Electrons in Crystal 
Lattices. Proceedings of the Royal Society of London. Series A, Vol. 130, No. 814 
(Feb. 3, 1931), pp. 499-513. 

[14] E. Landau, Uber die Einteilung der positiven ganzen Zahlen in vier Klassen nach der 
Mindeszahl der zu ihrer additiven Zusammensetzung erforderlichen Quadrate. Arch. 
Math. Phys. 13, 305-312, 1908. 

[15] J. Marklof and Z. Rudnick, Almost all eigenfunctions of a rational polygon are uni- 
formly distributed. J. of Spectral Theory. 2 (2012), 107-113. 

[16] F. Oravecz, Z. Rudnick and I. Wigman, The Leray measure of nodal sets for random 
eigenfunctions on the torus, Annales de ITnstitut Fourier 58 (1), (2008), 299- 335. 

[17] S. Rahav and S. Fishman. Spectral statistics of rectangular billiards with localized 
perturbations. Nonlinearity 15 (2002), no. 5, 1541-1594. 

[18] M. Rccd and B. Simon. Methods of Modern Mathematical Physics, Vol. 2: Fourier 
Analysis and Self-Adjointness. Academic Press, London, 1975. 

[19] P. Seba. Wave Chaos in Singular Quantum Billiard. Phys. Rev. Letters, Vol. 64, No. 
16, 1855-58, 1990. 

[20] T. Sliigcliara, Conditions for the appearance of wave chaos in quantum singular sys- 
tems with a pomtlike scatterer. Phys. Rev. E 50, 4357-4370 (1994). 

[21] A. Snirel'man Ergodic properties of eigenfunctions. Uspekhi Mat. Nauk 29, no. 6 
(180), 181-182, 1974 

[22] J. G . van der Corput, Neue zahlentheoretische Abschdtzungen. Math. Ann. 89 (1923), 
no. 3-4, 215-254. 

[23] S. Zclditch. Uniform distribution of eigenfunctions on compact hyperbolic surfaces. 

Duke Math. J. 55, no. 4, 919-941, 1987. 
[24] S. Zelditch and M. Zworski, Ergodicity of eigenfunctions for ergodic billiards. Comm. 

Math. Phys. 175, 673-682, (1996). 

Raymond and Beverly Sackler School of Mathematical Sciences, Tel Aviv 
University, Tel Aviv 69978, Israel 

E-mail address: rudnickOpost. tau.ac.il 

Raymond and Beverly Sackler School of Mathematical Sciences, Tel Aviv 
University, Tel Aviv 69978, Israel 
E-mail address: henrik@post.tau.ac.il 



